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fo3iowlz»g  ej^rs  have  T&eea  noted  in  tlie  derivations  of  tlie  fomal&s 
ap^e^j^  in  Part  11  of  EADC-TDR-6Jj-33^>  "Design  of  Metal  Space  Frame  Radomes^" 
dated  Rovembef  Since  tliey  occdr  in  tile  derivations  tlielr  influence  is 

refiecr^d  ttrpu^ut  Part  II  and  tiielr  correction  will  require  >STibstantial 
rewriting#-  The  revision  will  hopefully  te  available  by  May  1965» 

R 

Page  75>  Equation  (3).  This  equation,  b77R^g:=  l/2vd^  L#  should 

1 

read  iiTTR^ktg-s  l/2wd  ^  to  account  for  truncation  of 
the  ^hera. 

Page  75,  Equation  (5)«  This  eqj^atlon  tg=;  if  d  should  read 


d. 


Page  82,  Equation  (3l)»  term  P  should  be  twice  the  density  of  the 
actual  mterlal  raed  in  the  r^dome  members  since  tg  is  representative  of  half 
the  actual  volume  of  material  in  the  members. 


Page  86,  Equation  The  R^  in  the  denominator  of  the  second  term 

of  this  equation  should  be  // 

Page  86,  Equation  (^i-l).  This  equation  is  based  upon  an  Euler  backing 
stress  for  a  fixed  ended  column.  Since  you  cannot  achieve  complete  fixity  at 
the  ends  of  the  radome  members  this  is  considered  too  liberal  an  assumpticn. 
The  more  conseanrative  approach  would  te  to  assume  the  members  as  pin  ended. 
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T?valuation 


This  contractual  effort  resulted  in  the  develcroment  of  two  classes  of 
theoretical  tools  heretofor  not  available  to  the  radesne  designer  - 

1.  Simplified  formulas,  charts  and  graphs  for  coaputing  the  insertion 
loss,  side  lobe  perturbation,  and  boresi^t  shift  caused  by  a  metal  space 
frame  radone.  Ifliereas  these  co^utations  foimerly  required  extensive  time 
and  can5)uter  ser\’’ices,  they  now  may  be  performed  in  a  timely  fashion  by 
ordinary  desk  cdcputatibnal  techniques.  This  allows  us  to  make  aimlytical 
checks  oh  radome  performance  which  until  this  time,  have  not  been  possible 
bn  a  practical  working  basis. 

2.  An  integrated  electrcmiagnetic  -  structural  design  procedure  for 
metal  space  frame  radanes.  The  design  is  conducted  within  simulteiheous 
structural  and  electromagnetic  constraints  resulting  in  a  presc'»'ibed 
electrical  transmission  loss  coupled  with  a  prescribed  level  of  structural 
integrity  at  minimum  practical  wei^t.  Previbus  radome  design  procedure 
has  been  ah  iterative  process  of  structural  design  v.s.  electrcanagnetic 
checks  mtil  ah  acceptable  design  was  achieved. 

The  twb  techniques  described  above  represent  a  large  step  in  a  long  range 
program  to  convert  ground  radome  design  to  a  practical  engineering  procedure 
which  is  not  dependent  upon  extensive  ccmputer  programs  and  electromagnetic 
testing. 

This  information  will  be  made  available,  through  distribution  of  the  final 
report j to  those  agencies  and  firms  involved  in  the  design  and/or  application 
of  rigid  grbund  radomes. 
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ABSTRACT 


PART  1  -  Design  graphs,  formulas  and  procedures  for  the  design  and  evalu¬ 
ation  of  the  electrical  performance  of  metal  space  frame  radomes  is  pre¬ 
sented  covering  specifically  transmission  loss«  boresight  error,  and  side- 
lobe  variation.  The  relative  merits  of  metal  and  dielectric  elements  is 
analyzed. 

PART  II  -  The  study  consists  of  developing  a  straight  forward  simple  pro¬ 
cedure  for  the  optimum  sizing  of  solid  beam  elenients  for  large  space  frame 
radomes.  Equations  were  derived  e^ressing  the  relationship  between  the 
element  ditnensions  necessary  to  satisfy  a  prescribed  electrical  tr^sinissioh 
loss  and  the  structural  integrity  of  the  space  frame.  These  equations,  con? 
sidered  as  design  constraints,  were  utilized  in  the  selection  of  the  element 
size  which  satisfied  a  defined  optinium  criterion.  Two  bptirhization  criteria 
were  adopted:  first,  the  minimization  of  the  transmission  loss  subject  to  a 
constraint  of  structural  integrity  and,  second,  the  reduction  of  the  total  struc¬ 
tural  weight  of  the  elements  which  simultaneously  satisfy  the  transmission 
loss  and  structural  constraint. 

A  sample  design  computation  was  conducted  as  an  example  of  the 
procedure.  The  second  optimization  criterion  was  applied  to  a  150  foot  diam¬ 
eter  radome  which  was  required  to  sustain  150  mph  winds.  A  procedure  for 
implementing  the  first  optimization  criterion  employing  the  method  of 
Lagrangian  multipliers  was  presented  resulting  in  four  equations  to  be  solved 
simultaneously. 
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PART  I 


SYMBOLS 

(also  used  in  Appendix) 


G  -  relative  power  loss  due  to  space  frame  radbme 

p  -  blocking  area  ratio  of  space  frame  radons  >  i.  e.  ratio  of  projected 
area  of  space  frame  to  area  of  reflector 

g  -  the  induced  field  ratio  of  a  scatterer,  i.  e. ,  the  field  strength  in 
the  projected  blocking  area  of  the  scatterer  which  would  produce 
the  actual  scattered  far  field  produced  by  the  scatterer,  divided 
by  the  incident  field  strength 

W  -  width  of  scattering  element 

d  -  depth  of  scattering  element 

L  -  average  length  of  scattering  element 

g^  ,  g^  *■  values  of  g  when  the  polarization  of  the  incident  field  is 

respectively  longitudinal: and  transverse  toihe  length  of  the  element 

X.  -  free  space  wavelength 

f  -  frequency  of  electromagnetic  radiation 

c  -  parameter  depending  on  curvature  of  radome 

D  -  hub  diameter  also,  in  context,  the  dish  aperture  diameter 

P  blocking  area  ratio  of  the  hubs 

c,  £  -  dielectric  constant  relative  to  free  space  value  and  thickness, 
respectively,  of  skins 

r,r^  -  polar  coordinates  in  the  antenna  aperture 
^  (r,  /()  -  pbAse  error  at  (r,A| ) 

0,  O'  -  boresight  shift  in  radians  due  to  a  phase  error 
Orms  ~  bdresight  shift 
a  -  radius  .of  dish  aperture 

f(r)  -  illiunination  of  circularly  symmetric  aperture 
/^g«  ,  ^L,  ^d,  -  variations  in  these  quantities 

[rt  j  ~  I  numb^<r^|  elements,  hubs,  or  windows  having  a  specified 

tolerance  or  design  variation 
(g/tg,)/2 

*^1’  *^1*’  ^2’  ^3’  ^4’  ^5’  ~  ^^ctions  of  aperture  illumination 

defined  in  context 


A  -  area  of  a  part  of  the  radome 
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SYMBOLS  (cont'd) 


25.  R  -  "radius  of  radome 

26.  0,t  -  defined  in  Figure  38 

27.  0^  -  angle  at  which  grating  lobe  occurs 

28.  tp  -  defined  in  eq(34) 

29.  P  -  defined  in  eq(37) 

30.  A  ,A  -  defined  in  eq(36) 

e  o 

31.  subscript  m  -  metal 

32.  '*  d  -  dielectric 

33.  E  -  Young's  modulus 

34.  S',  a,  <r  ,  (T.  ~  quantities  defined  in  |3j 

m  d 
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PART  II 
SYMBOLS 


1.  A  -  area 

2.  C  -  constant 

3.  D  -  constant 

4.  d  -  element  depth 

5.  £  -  modulus  of  elasticity 

6.  F  -  total  force 

7.  f  -  unit  force 

8.  H  -  constant 

9.  I  -  moment  of  inertia 

10.  k  -  surface  area  coefficient 

11.  K  -  Tslen-Von  Karman  coefficient 

12.  L  -  element  length 

13.  M  -  moment  or  number  of  panels 

14.  N  -  stress  resultant  or  nvunber  of  elements 

15.  Pq  -  pressure  at  wind  stagnation  point 

16. ~  q  -  distributed  transverse  load 

17.  Q  -  equivalent  shell  weight  density 

18.  R  -  radius  of  radome 

19.  S  -  effective  spacing  of  elements 

20.  t  -  thickness 

21.  u  -  circumferential  displacement 

22.  V  -  velocity  of  wind 

23.  V  -  meridional  displacement 

24.  W  -  weight  of  elements 

25.  w  -  element  width 

26.  X  -  coordinate 

27.  y  -  coordinate 

28.  z  -  coordinate 

29.  Cl  -  meridional  angle  from  base 

30.  y?  -  exponent  coefficient 

31.  0  -  circumferential  shell  coordinate 
32  X  -  Lagrangian  Multiplier 

33.  p  -  Poisson's  ratio 

34.  p  -  weight  density 

35.  normal  stress 

36.  ^  -  meridional  shell  coordinate 


37. 

w  -  normal  displacement 

38. 

subscript 

1  -  boundary  condition  one 

39. 

If 

2  -  boundary  condition  two 

40. 

If 

R  -  radome 

41. 

M 

r  -  reinforcement 

42. 

tl 

s  -  shell 
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SYMBOLS  (cont’d) 


subscript  e  -  external 
'•  i  -  internal 

"  cr  -  critical 

"  yp  ~  yield  point 

"  b  -  base  of  radome 

"  ^  -  direction  of  meridian  coordinate 

"  0  -  direction  of  circumferential  coordinate 

'•  ^0  -  direction  of  ^  in  plane  normal  to  0 

”  n  -  direction  normal  to  surface  of  radome 


DESIGN  OF  METAL  SPACE  FRAME  RADOMES 
PART  i  -  ELECTRICAL  DESIGN 


INTRODUCTION 

"Electrical  Design  of  Metal  Space  Frame  Radomes-  '  £l]  by  Alan  F.  Kay, 
was  prepared  and  submitted  for  publication  on  August  20, 1963^  A  revised  ver¬ 
sion  dated  March  23,  1964  has  been  accepted  for  publication  in  the  Transactions 
of  the  IEEE,  Professional  Technical  Group  on  Antennas  and  Propagation.  This 
paper  contains  the  fundamental  theory  on  which  most  of  this  present  report  is 
based.  The  present  report  gives  design  graphs,  formulas,  and  procedures  for 
specific  radomes  in  the  .  1  to  lOgc/s.  band,  as  well  as  considerably  more  design 
data  than  included  in  £1] .  This  data  pertains  to  the  three  most  critical  electri¬ 
cal  properties  of  the  radome:  loss,  boresight  error,  and  sidelobe  increase. 

A  section  on  the  relative  merits  of  metal  and  dielectric  ribs  is  included. 

LOSS  DUE  TO  SPACE  FRAME  ELEMENTS 

Equation(22)  of  [ij  is  the  fundamental  approximate  formula  for  the  loss 
due  to  the  space  frame  itself.  In  the  Appendix,  p.  120,  a  brief  derivation  is  given 
of  this  expression: 

(1)  G=|l  +  /’gP 

where  G  is  the  relative  power  loss  as  measured  at  (or  near)  the  peak  of  the 
antenna  beam  due  to  the  presence  of  the  radome  space  frame,  p  is  the  block¬ 
ing  area  ratio  of  the  space  frame  and  g  is  the  induced  field  ratio  of  the  space 


1 


frame  elements.  £qu«ition(l)  neglects  the  loss  due  to  hubs  and  dielectric  skin. 
We  give  here  simple  formulas  and  graphs  by  which  G  may  be  computed  in 
most  cases. 

Consider  the  space  frame  to  be  comprised  of  triangles  of  elements  of 
average  length  L  from  hub  center  to  hub  center  and  of  width  W.  Consider 
any  one  triangle  as  extending  to  the  midline  of  each  of  the  three  elements  com¬ 
prising  it.  If  the  triangles  were  exactly  equilateral  then  P  would  equal 


W 

If  the  triangles  were  right  isosceles  then  ^  =  3. 56  .  f  generally  in¬ 

creases  the  further  from  equilateral  the  triangles  depart  or  the  greater  the 
standard  deviation  of  the  element  lengths  from  the  average  value  L.  For  a 
well  designed  radome  geometry  the  element  lengths  should  all  be  close  to  the 
average  L  and  the  triangles  all  reasonably  close  to  equilateral.  A  reasonable 
approximate  value  of  P  has  been  found  to  be 

(2) 


This  value  is  also  approximately  corrected  for  a  quadrilateral  geometry. 

The  value  of  g  depends  on  the  element  croas^ection,  the  average  inci¬ 
dence  angle  0  of  the  radiation  from  the  dish  to  the  curved  radome,  and  the 
polarization.  When  the  projection  of  the  elements  into  the  aperture  plane  are 
all  parallel  to  E,  the  appropriate  value  of  g  is  g^  .  When  the  projection 


2 


of  the  elements  into  the  aperture  plane  are  all  transverse  to  E,  the  appropriate 
value  of  g  is  g^.  In  the  case  of  a  triangular  geometry  with  L.  <  <  D, 
where  D  is  the  dish  diameter,  the  appropriate  value  of  g  is 


go 


8av  = 


Values  of  g  for  various  rectangular  cross  sections,  W  x  d,  and 
frequencies  are  shown  in  Figures  1  to  3  for  0  =  0,  These  values  were  ob¬ 
tained  by  the  computer  program  described  in  At  sufficiently  high  fre¬ 
quencies  f,  ^{g  )  approaches  -1  and  ^*/(g  )  approaches  zero.  How- 

2iV  SV 

ever,  the  manner  of  approach  is  such  that  each  of  the  curves  d/W  =  constant 
spirals  an  infinite  number  of  times  before  reaching  the  terminal  point  (-1, 0). 

If  the  graphs  were  extended  into  the  regime  where  this  spiraling  takes  place, 
t  hey  would  be  difficult  to  read  and  interpolation  either  in  f  or  d/W  would 
be  virtually  impossible.  It  is  also  in  this  regime  where  shaping  of  the  element^^^ 
can  effectively  reduce  loss  and  consideration  of  purely  rectangular  elements  be¬ 
comes  less  significant.  Figures  4  to  6  are  transparent  overlays  for  three  values 
of  /'  which  give  G  directly  in  db,  according  to  equation(l)  if  any  of  these  are 
placed  over  one  of  the  Figures  1-3  with  the  axes  coincident.  In  any  case  where 
values  of  /?,  W,  and  d  imply  that  L  <  X,  then  multiple  scattering  becomes 
important and  equation(l)  is  inaccurate.  These  cases  occur  at  low  f  when  the 
apparent  values  of  G  obtained  from  (1)  seem  to  increase  above  unity  (negative  db) 
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with  decreasing  f.  Figures  7-10  show  the  loss  of  a  particular  radome  versus 
frequency  for  various  average  incidence  angles  0.  0  =  0  corresponds  to  a 

flat  radome.  Typically  Q  is  10  to  20*^  for  spherical  radomes.  Each  figure 
is  for  different  element  widths  and  depths  as  indicated.  A  similar  plot  was 
shown  in  Figure  10  of  Q]  for  still  a  different  element  including  comparison 
with  measured  values.  Figure  11  is  a  photograph  of  measurements  of  inser¬ 
tion  loss  of  a  partial  radome.  Figure  12  shows  measured  and  theoretical  loss 
of  two  such  radom,es,  which  were  designed  for  optimum  electrical  and  struc¬ 
tural  performance  as  described  in  Part  II,  p.  100. 

An  empirical  formula  for  G  has  been  obtained  by  observation  of  the 
general  behavior  in  Figures  1  to  3  and  other  similar  data,  as  follows! 

(4)  101ogjQG=  i  [39.5W  +  3.5d+.l47^  \ 

where  c  is  a  parameter  depending  only  on  the  curvature  of  the  radome  over 
the  projected  area  of  the  antenna,  or  approximately  on  the  average  incidence 
angle  0.  If  0  =  12°.  c  =  4.46,  If  0  =  25°,  c  =  11.8.  Equation(4)  is  reasonably 
valid  provided 

W  <  2\,  W<d<10W,  1 10  log^^Gj  <  3  db. 

This  validity  is  shown  in  Figures  13-29  where  equaticn(4)  is  compared 
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with  equation(l}.  From  eq(4)  we  observe  that  depth  d  appearing  in  the  second 
term,  is  about  one -tenth  the  importance  in  producing  loss  of  fte width  w  appear¬ 
ing  in  the  first  term.  The  third  term  expresses  the  frequency  dependence  and 
shows  that  the  loss  of  a  given  radame  tends  to  increase  with  increasing  X.  at  a 
rate  proportional  to  d/w.  The  last  term  expresses  the  effect  of  curvature  of 
the  radome.  The  latter  effect  becomes  negligible  as  the  element  becomes  square 
(d  — >  w)  or  at  sufficiently  low  frequencies  or  low  average  incidence  angles.  In 
any  case  the  first  term  usually  dominates  this  expression.  InFigures  13-29. 
the  approximation  (4)  to  the  more  exact  formula  (1)  is  probably  more  accurate 
for  the  cases  where  the  two  disagree  considerably,  since  this  disagreement 
occurs  when  mutual  coupling  is  important  and  the  loss  indicated  by  (1)  is  too 
small. 

LOSS  DUE  TO  HUBS 

Let  D  be  the  hub  diameter  and  L  the  average  element  length  from 
hub  center  to  hub  center.  The  analogous  formula  to  (1)  including  the  effect  of 
hubs  is 

(5) 


where  is  the  blocking  area  ratio  of  the  hubs,  assuming  optical  blocking 

of  the  hubs.  In  a  triangular  or  quadrilateral  geometry  is  approximately 


(6) 


71  py 


D.2 


p  =  ^  Q/;:) 
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Hence  (5)  can  be  written  as 


(7) 


G  =  1  + 


3.5W 


£ 


L  +  D 


D  <  <  L. 


With  reasonable  structural  design  the  effect  of  the  hubs  is  usually  a  small 
correction  to  the  effect  of  the  elements. 

LOSS  DUE  TO  DIELECTRIC  SKINS 

In  a  space  frame  radome,  at  least  up  to  lOgc,  the  skins  will  be  thin 
compared  to  a  wavelength  and  the  loss  will  be  entirely  due  to  reflection  for 
practical  purposes  if  low  loss  dielectric  windows  are  used.  The  reflection 
loss  in  db  is 


(8) 


10 


,,  ^2  .  Z,Z-nS>sr €  4 

(l-€)  sin  ( — - ) 


where  ^  is  the  skin  thickness  and  €  is  the  dielectric  constant.  Equation(8) 
applies  to  normal  incidence  but  is  a  good  approximation  for  spherical  radomes. 
For  thin  radomes  a  sufficient  approximation  is 


(9) 


10  1ogjQ(lM«-l)^^)^) 


For  radomes  below  2gc  the  loss  expressed  by  (9)  is  usually  entirely 
negligible.  Above  2gc  it  rises  rapidly.  Sufficiently  thin  skinned  radomes  can 
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only  be  realized  above  lOgc  with  considerable  care  in  structural  design. 
BORESIGHT  ERROR 

Equation  (26)  of  {]lj  is  a  fundamental  approximate  formula  for  boresight 
error  due  to  a  space  frame  radome  or,  if  properly  interpreted,  for  any  source 
of  aperture  phase  errors.  This  equation  is  as  follows: 

a  2:r 

.  J.rf(r)^f(r.n)  cos/J^dr 
(10)  O'  =  - 

ir^f(r)dr 

Here  0*  is  the  boresight  shift  in  radians  in  the  plane  0  where  f(r)cos/)^ 
is  an  approximation  to  a  monopuXse  difference  pattern,  (r,/|^)  are  polar  co¬ 
ordinates  in  the  aperture  of  radius  a,  and  X.  is  the  wavelength,  is 

the  phase  error  at  (r,^).  In  this  section  we  shall  use  this  equation  to  determine 
the  boresight  error  due  to  many  commonly  occurring  factors.  Equation  (10)  and 
the  results  of  this  section  apply  to  monopulse  systems.  However,  they  are  rea¬ 
sonably  accurate  for  conical  scan  or  peak  shift  determinations  if  f(r)  is  con¬ 
sidered  as  the  aperture  illumination  of  a  sum  or  even  symmetry  mode. 

The  types  of  errors  considered  are  indicated  schematically  in  Figures 
30-36  as  variations  from  the  small  boresight  shift  of  a  basic  radome  or  from 
no  boresight  shift  at  all. 

The  following  notation  is  used: 

X  =  wavelength 

a  =  radius  of  antenna  aperture  =  D/2 
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f(r),  0  <  r  <  a  =  illumination  of  aperture  (circularly  symmetric) 

W,  h,  d  =  nominal  element  width,  average  length,  and  depth  respectively 
^W,  ^d  =  tolerance  or  design  variation  in  W,  L,  or  d,  respectively, 

for  elements  specified 

=  blocking  area  ratio=T3.  5  W/L  for  triangular  or  quadrilateral  division 

of  sphere 

.^^g)  ±  imaginary  part  of  the  induced  field  ratio,  IFR.  If  polarization 
is  random  or  averages  to  45°  over  particular  elements  causing  boresight  shift, 
then 

-  'Sav 


is  the  algebraic  average  of  the  longitudinal  and  transverse  polarization  IFR's, 
l£  not,  applicable  g^  or  g^  or  weighted  average  must  be  used. 

~  tolerance  or  design  variation  in  J^ig)  of  elements  specified. 

A  =  total  area  of  radome  where  tolerance  or  design  variation  may  occur. 

of  elements,  hubs,  or  windows  having  a  specified 
tolerance  or  design  variation?  0  <  p  <  1. 


fraction 

In! 

number 

^  =  nominal  dielectric  constant  and  skin  thickness  of  windows 


/\c ,  =  variation  in  c  or  i' 


All  boresight  errors  are  in  radians. 
Note  that 


A(WJ!(g)L)  -  wX(g)^L  +  WL^,(g)  t-  L^.ig)^W  . 
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We  now  consider  six  different  sources  of  boresight  error: 

1)  A  small  hiifhber  rii  n  =  l,  2i  i,,  of  extra  elements  are  located  in  one  local 
region  of  the  radome,  small  compared  to  the  antehha'sdiameter  D.  The  worst 
error  occurs  when  the  extra  elements  are  at  the  peak  of  &e  difference  pattern, 
when  the  beam  shifts  away  from  this  location  by  ^  angle 

(11)  O'  =XhLw4^g)cQ 

where  g  is  the  IFR  value  of  the  extra  elements  and  c-  depends,  only  dh  the 
illumination  taper 

max  (l(r)^  (.  27, fa^,  -  lOdb  taper 

Fig.  30) 

2}  Certain  elements  of  the  space  frame  whic:h  are  randomly  distributed  in  the 
space  frame  and  constitute  a  fraction  p.  of  the.  total  number  of  elements,  have 
a  different  width  AW  and/or  depth  d^  =  +  ^d  from  the  remaining 

elements  which  are  assumed  to  be  such  that  if  =  0  then  there  is  no 

boresight  error.  RMS  boresight  shift  is  then 


(12) 


*^0  ■ 


0<r<a 


2ir^5r^f{r)dr 


.37 /a  ,  20db  taper 


(see 


(13) 


O' 

rms 


)VV  Cj 


triangulated 

■ 

radome 

V 
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where 


(14) 


^1  = 


i  ^  2 

j  3.  5Tr|  rf  (r)dr 
2if^^r^£(r)dr 


,42/^  lOdb  taper 


Sd/a  ,  ZOdb  taper 


K  the  radome  is  not  triangulated  by  the  space  frame,  but  has  some  other  arrange 
merit  of  elements,  (12)  applies  wi&.  p.  replaced  by 


(15) 


3.5W 


(see  Fig. 31 )  , 


3)  Given; tolerances  and  ^d  in  dirhensions  of  all  elements 

a)  Largest  possible  boresight  error  occurs  when  +^W ,  +^d  occurs 
on  one  side  of  dish;  -^W,  -^d  occurs  on  other  (Figure  32a).  Beam  shifts 

towards  the  least  mstal  area  by  an  aihount 

)^A(W«Sig))c2 

(16)  9  = - r- - 


where 

(17) 


a 

7^rf(r)dr 

•iT^ijr^f(r)dr 


I  1.19/a,  lOdb taper 
\  1.30/a,  20db  taper 


b)  Variation  of  ,  Ad  is  limited  to  an  area  A  whose  maximum 
diameter  is  small  compared  to  D.  Boresight  shift  is 
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where 


XA^wijig))c 


3.  5~m^  |^^  f(r)  ‘| 
_ Q<^a_.  J 

2-  2 

Zf  ^r  f(r)dr 


.  947a  j  i  Odb  .toper 


a;  1,  36/a^i 


ZOdb-  taper 


(aee  Fig.  32b.). 


c):  Variations  in  dimensions  of  elements  of  a  space  frame  are  random  with 
and  ^d  being  rms  values  of  the  variations  in;  W  and  d.  Then  rrns  value 
of  boresight  error  is 


^rms  "  (see  Fig.  32c) 


I  4)  Rms  boresight  error  for  a  random  space  frame  geometry  (most  probable 

I 

I  boresight  error  if  a  number  of  elements  are  randomly  distributed  over  sphere) 


(21a) 


0  =  XW.5l,(g)LVnc  > 

rms  ^  ®  1 


‘=l’  = 


1  Jjf^(r)r  dr 
2^2 

Zif  a^r  f(r)dr 


I .13/a^, 


lOdb  taper 


.  15/a  ,  20db  toper 


where  h  =  ira  p  J'WL  is  the  number  of  elements  seen  by.  the  antenna.  For  a 
triangular  or  quadrilateral  geometry 


P- 


and  then 


(21b) 


0  . 
rms 


(see  Fig.  33); 


Mutual  coupling  between  elements  is  neglected  in  {21a); 

5)  A  small  number  n,  n  =‘  1 , 2, , .  *  of  dielectric  windows  (each  of  average 
area:  4L  )  with  variations  and  in  dielectric  Constant  and  wall 
thickness  are  all  located  in  an  area  of  the  radoine  whose  maximum  diameter 
is  small  compared  to  a.  The  maximum  boresight  error  occurs  when  these 
windows  are  at  a  difference  pattern  peak,  and  is  given  by 


(22) 


0  =  nL  c 


(V€-l)/^^  + 


(23) 


^^4^  a  2 


.  4  max  |f(r)|| 
0<r<a 


■>fSr  f(r)dr 


j  . 68/a^,  lOdb  taper 

3 

.94/a  ,  20db taper 


V 


(see  Fig.34) . 
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b)  Variation  is  random  and. windows  are  distributed  uniformly  over 


surface  of  radome  with  and  ^  bein;»  rms  values  of  the  variation  in  c 
and  ^  .  Then  rms  value  of  boresight  error  is 


(27) 


rms 


('sTc-l)^  + 


Lcg  (see  Fig  361^. 


fom 


Figure  37  is  a  plot  of  ■T^\^g  )  for  various  d/X  and  W/X.  Frc 

\  SiV 

this  graph  the  preceding  formulas  may  be  evaluated  explicitly.  As  an  exaniple, 
we  consider  two  cases  -  the  first  a  gross  geometry  with  loose  tolerances,  the 
second  a  finer  geometry  with  tighter  tolerances. 
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The  resulting  values  of  boresight  shift  computed  from  eqs{14)  to  (27) 
are  shown  directly  on  the  figures,  30  thru  36  as  cases  1  and  2. 

EXPEL  :^ED  SIDELOBE  LEVELS  DUE  TO  THE  SPACE  FRAME 

As  pointed  out  in  pp.  14-15  of  £lj  the  expected  close-in  siddlobe  level 
due  to  the  presence  of  tlie  radome  increases  even  though  the  expected  average 
field  values  are  unchanged  by  the  radome.  A  rough  formula  for  the  expected 
increase  in  close-in  relative  sidelobe  level  relative  to  the  peak  of  the  pattern 
is 


(28) 


which  is  derived  in  the  Appendix, 

Values  of  average  sidelobe  level  increases  computed  by  (28)  are  usually 
small  compared  to  peak  sidelobe.  increases  which  may  even  be  an  order  of  mag¬ 
nitude  larger.  However  (28)  is  useful  in  showing  how  the  various  parameters 
effect  expected  sidelobe  increase. 

Another  sidelobe  phenomenon  often  occurs  with  space  frame  radomes, 
namely  at  some  angle(d)  far  from  the  main  lobe,  sidelobe(s)  may  arise  as  high 
or  higher  than  the  close -in  sidelobes.  There  are  called  the  diffraction  side- 
lobes,  in  anabgy  with  diffraction  grating  effect.  Since  to  the  extent  that  the 
elements  are  pa.-allel,  equi-spaced  and  coplanar,  i.  e.  a  diffraction  grating, 
these  lobes  will  be  higher  and  more  noticeable.  At  angles  out  to  0  <  0  <  \/2Li, 
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the  expected  field  values  with  the  radome  are  proportional  to  those  without, 
the  proportionality  constant  being  independent  of  0.  Except  for  the  snnall 
effect  considered  above  {or,  in  detail,  in  the  Appendix),  the  relative  sidelobe 
level  with  and  without  the  radome  in  the  range  0  <  0  <  X./2L  is  the  same.  If 
the  scattering  patterns  of  the  individual  elements  were  isotropic  outside  of 
this  range  and  if  the  scattered  fields  were  randomly,  or  even  better,  uniformly, 
distributed  in  phase,  the  energy  scattered  by  the  space  frame  would  be  scattered 
essentially  isotropically,  i.  e.  uniformly  distributed  in  all  space  and  the  far- 
out  sidelobe  levels  would  be  well  below  the  isotropic  level  of  the  total  radiation 
and  hence  usually  not  observable  except  on  an  average  basis,  by  sensitive  radio- 
metric  measurements.  However,  this  is  not  generally  the  case  and  diffraction 
lobes  do  arise. 

♦ 

The  IBM  7090  program  will  predict  these  lobes  accurately.  However, 
in  order  to  obtain  an  approximate  formula  for  rapid  calculation  of  their  levels 
and  to  see  how  and  to  what  extent  they  will  occur,  we  derive  here  an  approximate 
formula  for  them. 

Suppose  Li  is  the  average  length  of  an  element,  and  suppose  the  direction 
of  propagation  is  along  the  z-axis.  Suppose  an  element  is  at  a  distance  t  from 
the  axis.  The  phase  ^  as  observed  in  direction  0  of  the  field  incident  on  the 
element  is 

(29)  k(OB  -  AC)  =  Rnfl  -  (!-)"■  -  kR  co8(0  -  sin“^^) 

*  described  in 
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■(30) 


If  t  < 


(31) 


Referenced  to  the  axial  point  t  =  0,  this  becomes 


^  ^  =  kR(*fi-{-^^  -  1  +  cos  0  -  cos(0 


sin  (— )). 


<  R,  this  is  approximately 


2 

A  "  hR(-^sin  0 - -  kt  sin  0 


Now  if  the  radome  were  flat,  the  last  term  in  (31)  would  be  zero  and,  it  turns 


out,  (31)  would  be  exact 


% 


(32)  ^  ■  Aq  "  ® 

Hence  the  last  term  in  (31)  is  the  first  order  correction  term  for  phase 
variation  of  the  fields  scattered  by  the  elements  due  to  the  curvature  of 
the  radome. 

If  L  is  the  average  element  length,  the  average  element  center  to 
center  spacing  is  abou*  L,/2  for  a  triangulated  space  frame  (see  Figure  39). 


Figure  indicates  distance  between 
phase  centers  averages  about  L«/2. 

FIGURE  39 

To  a  first  approximation,  in  the  case  of  a  flat  radome,  from  eq(32)  a 
diffraction  grating  lobe  would  occur  at  an  angle  0  where 


(33) 


kL 

2 


sin  0 

g 


-  n. 


The  curvature  correction  term  in  (31)  implies  that  only  a  certain  number 
of  elements  near  the  axis  contribute  coherently  to  the  grating  lobe,  those  ele¬ 
ments,  in  fact  for  which  t  is  less  than  t^j  where 


(34) 


=  ir  . 
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Other  elements  for  which  t  exceeds  the  limit  set  by  (34)  are  essentially 
incoherent  or  have  random  phase.  From  another  point  of  view,  elements  with 
t  <  tp  are  in  the  first  Fresnel  zone  for  in-phase  contribution  to  the  scattered 
field  in  direction  0  =  0. 


g 


The  gain  G  (0  )  of  the  scattered  field  in  direction  0  =  0,  equals  I/t 
^  g  g 

times  the  gain  of  a  uniformly  illiominated  disc  of  radius  t^  multiplied  by  the 


fraction 


/l«av 


of  the  relative  effective  blocking  area  of  the  space  frame. 


multiplied  by  2,  and  finally  multiplied  by  the  periodicity  factor  P 


(35) 


2 

G  (C  )  =  8/Jlg  j  Pu  (if  t  <  <  R). 

s  g  >  '  av  0 


The  factor  2  is  an  approximation  to  the  fact  that  the  elements  near  the  center 

of  the  aperture  are  more  heavily  illuminated  than  an  average  element  in  the 

aperture  because  of  the  antenna  illumination  taper.  The  factor  1/it  is  due  to 

the  phase  error  in  the  first  Fresnel  zone  (which  has  a  max  of  ir,  by  definition). 

The  validity  of  the  use  of  this  factor  is  proved  in  ref  >  p.  15,  case  N  =  1. 

The  periodicity  factor  P  is  a  measure  of  the  parallelness  and  periodicity  of 

the  elements.  P  may  be  defined  as  follows!  Consider  the  Fresnel  zones  in 

the  projected  dish  aperture  as  observed  from  the  angle  0  .  These  are  parallel 

g 

strips  of  width  equal  to  the  average  periodicity  c£  the  space  frame  (approx  L/2). 
Let  the  total  element  area  in  the  even  and  odd  zones,  respectively,  be 


(36) 


A  LW  and  A  ^  N  LW 
e  —  e  0—0 
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where  N  and  N  are  the  numbers  of  elements  in  the  even  and  odd  zones 
e  o 

respectively.  Then 


(37) 


P  = 


A  -  A 
e  o 


A  +A 
e  o 


If  all  of  the  elements  are  parallel  (diffraction  grating)  then  P  =  1  (Figure  43). 

H  the  elements  are  in  a  regular  rectangular  array  then  P  =  1/Z  (Figure  41). 

If  the  elements  are  in  a  regxilar  triangular  array  then  P  =  1/3  (Figure  42). 

For  a  regular  hexagonal  array  in  the  orientation  shown  in  Figure  43 

.6. 

The  preceding  considerations  must  be  modified  if  the  elements  in  either 
the  even  or  odd  Fresnel  zones  have  a  strongly  preferred  direction  which  is 
not  at  45°  to  the  polarization,  and  if  jg^  }  is  substantially  larger  than  jg^J. 

For  example,  in  the  case  of  the  hexagonal  arrays  or  rectangular  arrays, 
if  jg^l  >  >  jg^l  and  the  polarization  is  parallel  to  the  Fresnel  zone  boundaries, 
then  P  might  be  close  to  unity. 

Assuming  a  50%  aperture  efficiency  with  the  radome  in  place,  the  gain 
of  th<"  antenna  is 


(38) 


1  0  2 
^0  "2^-1-^ 


Combining  (34),  (35)  and  (38)  gives  G  (0  )  relative  to  the  maximum  gain  of 

8 

the  antenna. 
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Triangular 
Qtamatry  with 
Fraaaal  Zonas 
Sha4ad  for  Diffrac¬ 
tion  Orating  Lobo 


Oaomatry 
with  Fraanal  Zones 
Shaded  for  Diffraction 
Orating  Lohe 


r/GU^ 


(39) 


G  (0  ) 

-J—i-  = 


ppuj"- 


kR<  c 


R  =  oo 


We  may  make  a  table  of  maximum  grating  sidelobe  level  as  observed  in  the 
IBM7090  calculations  and  tho.'je  predicted  by  (39)  for  v&rious  cases  where 
the  quantity 


(40) 


16XR 


ttD 


=  .045  =  13.5db 


Element 

IBM7090  Calc. 

Etj 

uaticn  (39i 

Geometry 

No. 

Table  1 
of  [l] 

Max  Diff  Liobe 
(db) 

P 

lOlogjo 

Max. 

Diff. 

Lobe(db^ 

1)  Regular  Snub- 
dodecahedron 

52 

25.4  ' 

.33 

1.1 

26.0 

2)  Regular  Snub- 
dodecahedron 

6 

22.6 

.33 

5.7 

24.  0 

3)  Regular  Snub- 
dodecahedron 

53 

22.2 

.  33 

5.4 

23.  7 

4)  Random  Icosahedral 

52 

34.2 

.  0058 

8.  3 

34.  2 

In  the  regular  geometry  f'ase,  the  value  of  P  =  ,  33  was  chosen  because 
this  geometry  is  close  to  a  regular  triangular  caoe.  The  agreement  in  the  above 
table  shows  that  this  value  of  P  is  reasonable.  The  dependence  on  element  shape 
checks  out  fairly  well.  In  the  random  case  the  value  of  P  was  chosen  to  make 
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the  two  calculations  of  diffraction  lobe  level  agree.  This  result  shows  that  the 
random  geometry  elements  are  about  five  times  more  uniformly  distributed 
in  angle  than  the  elements  of  the  regular  geometry 
COMPARISON  OF  DIELECTRIC  AND  METAL  STIFFENING  RIBS 

In  the  design  of  a  space  frame  radome,  there  is  always  a  compromise 
between  structural  strength  and  transmission  loss.  The  question  arises  as  to 
whether  one  does  better  in  this  tradeoff  using  a  dielectric  or  metal  space  frame. 

In  each  case  we  would  normally  use  the  highest  strength  materials  available  - 
steel  for  the  metallic  space  frame  and  Fiberglas  laminate  for  the  dielectric 
frame.  If  the  transverse  dimensions  of  the  elements  are  comparable  to  a  wave- 

r2i 

length,  data  exists  which  indicates  that  the  scattering  cross '  section  or  IFR 

of  metallic  and  dielectric  elements  are  roughly  the  same.  Since  the  metal 

enjoys  a  considerable  strength  advantage,  there  seems  little  question  that  in 

this  case  the  metal  is  superior.  Let  us  consider  the  case  where  the  maximum 

diameter  is  small  compared  to  a  wavelength.  In  this  case,  there  is  no  advantage 

in  shaping  the  metal  element,  so  that  we  may  as  well  consider  it  rectangular 

oi  width  W  and  depth  d  .  There  is  even  less  dependence  on  shape  for  scattering 
m  m 

jiy  the  dielectric  element  so  that  we  may  consider  it  also  a  rectangle  of  width 
and  depth  d^.  T’here  are  two  rough  criteria  for  equivalent  structural  strength 
of  the  metal  and  dielectric  elements!  equal  exial  stiffness  and  equal  radial  bend¬ 
ing  stiffness,  either  of  which  may  apparently  govern  under  different  circumstances. 
These  imply  rv<ispectively. 
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(41) 


E  W  d  =EWd^ 
m  m  m  d  d  d 


(42) 


E  W  d  E^W  d/ 
m  m  m  d  d  d 


where  E  /E  is  the  ratio  of  metal  to  dielectric  Young's  modulus,  which  in 
m  d 

our  case  is  about  15.  Let  us  use  (41).  We  have  then  equivalent  elements  satis¬ 
fying 


(43) 


W  d  =15W  d  =A 
d  d  mm  d 


where  A  ,  is  the  section  area  of  the  dielectric.  If  (42)  were  used  instead, 
d 

very  little  difference  would  be  found  in  the  general  results  which  follow. 

r3i 

We  may  use  the  results  of  Mei  and  Van  Bladel^  to  compare  the  two 

structurally  equivalent  elements  electrically.  At  transverse  polarization,  from 

Figures  5  and  6  and  eq(7)  of  [l]]  ,  we  observe  that  lor  both  metal  and  dielectric 

the  narrow  dimension  should  be  transverse  to  the  direction  of  propagation  for  best 

results,  i.  e.  W  <  d,  so  that  Figure  6  is  appropriate.  It  can  be  observed  that 

typically  a  is  about  twice  as  large  for  metal  as  for  dielectric  with  the  Fiberglas 

2 

dielectric  .onstant  (^^  4)  and  that  the  term  2S  in  eq(7)  is  about  equal  to  or  smaller 
2  2 

than  the  term  a  /e^  for  metal.  Accordingly  the  scattering  cross  section  ratio 

(T  !  ^  for  elements  with  the  same  dimensions  is  about  4  or  5.  But  bot). 
m  d  - 
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and  are  proportional  to  the  square  of  the  section  area  in  wavelengths, 

and  hence  for  structurally  equivalent  elements  satisfying  (43)  the  ratio  is 


/ , 


(44) 


m 

^d 


4  or  5 
15^ 


Thus  the  meta!  ir.  !>v,perior  electrically  by  about  50  to  1. 

In  the  case  of  parallel  polarization  the  situation  is  different.  The 

asymptotic  formula  (10)  of  QlJ  can  not  be  used  since  it  becomes  inaccurate 

in  the  range  of  interest  (indeed,  blows  up  when  the  perimeter  of  the  metal 

rectangle  is  X/Zn).  However,  from  [j4j  eqs(15)  and  (16),  we  may  relate  the 

IFR  at  longitudinal  polarization,  g.  ,  to  cT  as 

X  m 


(45) 


(T  =k|g„W  {' 
‘^m  m* 


k  = 


Ztt 


F rom  eq(  U)ol  [3j,  with  a  dielectric  constant  of  e  in  this  case. 


146) 


2  3  2 

(€-1) 


and  from  (43).  (45)  and  (46)  cT  <  if 

m  —  u 


(47) 


2  3  2  2 

,  (€-1)  k  225W  d 

,  I  1 2  mm 

k  g  w  <  - - 

m*  -  4 
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or 


15{«-l)kd 

(48)  |g^  1  1  - 2 - ^ 

In  Figure  44  g  is  plotted  for  various  d^  and  frequencies. 

Where  equality  holds  in  (48)  jgj  is  much  smaller  than  |g^  j  and  hence 

15(£-l)kd 

(49)  ±  - 4  • 

We  have  plotted  where  equality  in  (49)  holds  in  Figure  44  lor  the  case  c  =  4. 
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PART  II  -  STRUCTURAL  DESIGN  AND  OPTIMUM  ELECTRICAL- 

STRUCTURAL  DESIGN 


SUMMARY 

The  object  of  this  project  was  thi:  construction  of  a  systematic  pro¬ 
cedure,  amenable  to  simple  calculation,  for  the  optimization  of  large  space 
frame  radome  design.  Consideration  of  three  design  aspects  (electrical, 
structural,  and  mechanical)  formed  the  basis  of  the  optimization  procedure. 
Specific  attention  was  directed  to  radomes  of  140-160  feet  in  diameter  used 
with  RF  systems  radiating  within  a  200-500  me  band. 

The  project  study  involved  three  phases!  (1)  the  application  of  an 
electrical  performance  evaluation  technique  to  predict  radome  transmission 
loss,  (2)  a  survey  of  the  radome  literature  far  structural  analysis  and  mechani¬ 
cal  design  concepts,  and  (3)  the  consolidation  of  pertinent  design  considerations 
into  a  design  procedure  for  an  optimum  radome  within  practical  limitations. 

The  evaluation  of  the  electrical  performance  of  a  space  frame  radome 
involved  the  prediction  of  the  transmission  loss  due  to  scattering  effects  caused 
by  induced  currents  on  the  elements  of  the  blocking  structure.  The  theory  re¬ 
lates  the  gain  of  the  radome  to  the  cross  sectional  shape  and  length  of  the  space 
frame  elements.  A  specification  of  the  maximum  allowable  reduction  in  gain 
fixes  a  holonomic  relationship  between  the  width,  depth  and  length  of  the  elements. 
This  relationship  was  defined  as  the  electrical  constraint  on  the  radome  element 
design. 

The  accumulation  and  survey  of  the  pertinent  literature  (see  Bibliography) 
was  a  prerequisite  to  the  formulation  of  a  suitable  structural  constraint  on  the 
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element  design.  In  the  literature  is  provided  a  complete  structural  design 
guide  for  the  sizing  of  the  radome  elements  which  lacks  the  inclusion  of 
any  electricail  constraint.  Therefore,  it  is  not  possible  to  quantitatively  dis¬ 
criminate  the  degree  of  compromise  inherent  in  a  particular  design  between 
electrical  perforniance  and  structural  integrity.  This  is  not  a  satisfactory 
design  condition  if  optimization  is  the  objective.  It  was  concluded  that  the 
simultaneous  imposition  of  an  electrical  and  structural  constraint  on  the 
analysis  of  the  elements  would  constitute  the  design  approach. 

In  addition  to  the  literature  survey,  consultation  with  several  radome 
developers  and  manufacturers  was  conducted.  The  consultations  included 
personnel  of  Lincoln  Laboratory,  Bell  Telephone  Laboratory,  Goodyear  Air¬ 
craft,  North  American  Aviation,  and  M.  I.  T. 

Other  efforts  in  space  frame  radome  development  are  being  directed 
toward  hardening  against  nuclear  blast  induced  dynamic  loads.  This  concept 
is  not  considered  in  the  subject  study  which  directs  attention  to  static  load 
conditions  consistent  with  normal  environments.  The  significance  of  the  har¬ 
dened  radome  development  to  this  study  exists  in  the  fact  that  substantiation  of 
the  general  shell  analogy  concept  of  analysis  has  continuously  precipitated.  If 
advances  in  fundamental  design  procedures  exist,  they  are  apparently  numerical 
in  nature.  The  STRESS(MIT)  and  STAIR  (LL)  routines  are  examples  of  digital 
computer  programs  for  analyzing  space  frame  structures  of  high  order.  The 
matrix  theory  of  structural  analysis  provides  a  powerful  computational  tool  but 
is  limited  by  computer  capacity  and  as  yet  programs  have  not  been  composed 
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which  display  general  design  capability.  Since  the  project  specifications 
require  a  design  procedure  propitious  to  hand  or  simple  machine  calculation, 
the  equivalent  shell  analysis  was  adopted  for  the  optimization  procedure  rather 
than  alternatively  pu~suing  advanced  computer  program  development. 

The  formulation  of  a  structural  constraint  consisted  of  comparing  the 
criteria  of  strength  and  stability  of  the  space  frame  which  assures  the  integri¬ 
ty  of  the  structure.  The  criterion  which  imposes  the  most  severe  design  condi¬ 
tions  was  selected  as  the  structural  constraint. 

The  concept  of  optimum  radome  design  requires  the  definition  of  a 
criterion  by  which  comparison  of  selected  designs  can  oe  made.  There  are 
many  criteria  relative  to  which  good  design  principles  can  be  directed.  Two 
such  examples  are!  (1)  the  minimization  of  transmission  loss  subject  to  a 
structural  constraint,  and  (2)  reduction  of  structural  weight  in  conjunction 
with  an  electrical  and  structural  constraint.  Examples  of  these  two  approaches 
are  developed  later  in  this  report.  The  problem  of  producing  the  best  radome 
design  within  practical  limitations  cannot  be  solved  by  a  single  direct  analytic 
approach.  It  is  first  necessary  to  select  the  element  size  subject  to  an  optimi¬ 
zation  criterion  and  second  r.o  generate  sound  mechanical  design.  This  study 
is  predominantly  concerned  with  the  optimum  sizing  of  elements  and  includes 
some  recommendations  and  goals  for  the  ensuing  mechanical  design  effort. 
ELECTRICAL  CONSTRAINT  ON  THE  ELEMENT 

For  the  purpose  of  optimum  design  of  a  radome  over  the  200-500  me 
band,  we  may  write  eq(4)  of  Part  I  as 
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10  log^^G  =  {39.5W  + 3.5d+  j  ,  db 

if  we  neglect  the  curvature  ‘■e-'rr  which  is  usually  small  at  lower  frequencies. 
This  is  tantamount  to  observing  that  at  the  low  frequency  end,  the  curves  of 
Figures  7  to  10  are  almost  independent  of  0.  With  this  approximation 
it  is  clear  that  the  loss  is  greatest  at  the  lowest  frequency  in  the  band.  Setting 
X  =  61"  (I94  me)  and  10  log^^G  =  .5  we  may  write  an  electrical  constraint,  such 
that  for  any  radome  satisfying 

L  >  79W  +  7d  +  -^ 

with  all  dimensions  in  inches,  the  loss  will  be  less  than  .  5db, 

STRUCTURAL  ANALYSIS  OF  SPACE  FRAME 
A.  External  Lv^ads 

It  is  recognized  that  the  most  significant  load  modalities  on  a  spherical 
radome  are  those  associated  with  kinetic  energy  of  the  incident  wind  and  the 
structural  weight.  Ice  and  snow  loads  are  of  minor  magnitude  and  experience 
indicates  that  snow  and  ice  build-up  seldom  occurs  on  large  radomes. 

The  analytic  representation  of  the  wind  pressure  distribution  adopted 
for  the  structural  analysis  Wcis  that  based  on  the  potential  flow  theory.  A  poly- 

r9i 

nomial  representation  of  the  pressure  distribution '■  ascertained  from  wind 
tunnel  data  is  mors  realistic  than  the  distribution  of  the  potential  flow  theory, 
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but,  since  little  deviation  between  the  two  formulations  exists  on  the  windward 
side  of  the  radome  and  significant  difference  exists  in  the  case  of  application  to 
the  shell  equations  the  potential  flow  representation  was  accepted.  The  fact 
that  buckling  stability  is  critical  only  on  the  compressive  or  windward  side  of 
the  radome  where  the  two  representations  agree  (see  Figure  1)  helps  justify 
the  choice.  Figure  1  illustrates  the  variation  in  magnitude  of  the  two  formula- 

[7] 

tions  for  typical  wind  tunnel  data. 

B.  Structural  Analysis 

The  equations  which  establish  the  conditions  for  stability  and  strength 
of  the  radome  were  developed  by  applying  classical  theory  of  elasticity  to  a 
spherical  isotropic  shell  of  structural  equivalence  to  the  radome  space  frame. 
This  technique  is  well  defined  in  the  literature  and  was  closely  followed  in 
the  development  of  a  structural  design  constraint. 

The  conversion  of  the  radome  to  an  equivalent  shell  requires  two  inde¬ 
pendent  equivalences,  elastic  and  stress.  The  elastic  equivalence  is  stated 
by  equating  the  extensional  and  flexural  stiffness  of  the  radome  to  that  of  the 
shell.  This  is  represented  as: 


The  stress  equivalence  is  required  due  to  the  inherent  difference  in 


elastic  behavior  of  a  space  frame  compared  to  a  shell.  The  stress  equivalence 


73 


relationship  equates  the  cross  sectional  area  of  an  element  to  a  ccmparitive 
cross  sectional  area  of  shell.  If  a  large  number  of  elements  make  up  the 
radome  structure  the  equivalence  condition  is  satisfied  when  the  volume  of 

r?! 

the  shell  is  one-half  the  volume  of  the  space  frame,  or  for  elements  of 
rectangular  cross  section. 


(3) 


4irR  t. 


5 


i 


Then,  since  the  element  resists  the  same  stress  as  some  equivalent  section 
of  shell 


(4) 


wd  =  Stg  . 


The  equations  expressing  the  equivalent  shell  thickness  and  modulus  of 
elasticity  are 


(5) 


tg  =  Vp  d 


(6) 


E  w 
R 

^^p  s 


In  order  to  apply  these  equations  the  effective  spacing  mutt  be  derived 
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I 

m 


in  terms  of  the  radome  parameters.  This  was  done  by  letting  M  equal  the 
number  of  panels  making  up  the  surface  of  the  radome  and  N  equal  to  the  total 
number  of  elements.  If  the  space  frame  is  considered  to  be  a  triangulated  struc- 
ture  typical  of  radomes  with  random  or  regular  geometries  then  3M  =  2N  and 
the  area  of  an  individual  panel  equals 


(7) 


A  = 
P 


6irR^k 

N 


where  k  is  a  constant  accounting  for  the  truncation  of  the  sphere.  The  value 
of  k  is  expressed  as  the  ratio  of  areas  of  a  truncated  to  complete  sphere  and 


is; 


(8) 


1  -  cos  ^ 


k  = 


.  i 


where  is  the  meridional  angle  to  the  radome  base.  The  average  panel 
area  was  also  expressed  in  terms  of  the  average  element  length  and  is 


(9) 


V3  2 
A  <  -^L 
p-  4 


depending  on  the  radome  geometry.  Since  the  larger  A^  is  the  more  conser¬ 
vative  the  design  must  be.  Equation(9)  was  taken  as  an  equality.  Equating 
equations  (7)  and  (9)  and  solving  for  the  total  length  of  elements  gives 

This  is  exact  only  for  a  non-truncated  radome.  For  small  truncated  ra¬ 
domes  it  is  an  approximation. 
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I? 


(10) 


NL  = 


24TTR^k 

•sTaL 


4’| 

II 

If 


The  effective  spacing,  S,  was  expressed  by  combining  the  equivalent 
stress  equations  (3)  and  (4)  such  that 


(H) 


S  = 


wd 


8i:R 


N 

—  1 
1 


and  substituting 


for  the  total  length  of  elements  gives 


(12) 


S  = 


8irR 


ZnTS  L 


(N-l)L  “  3(l-cos  ^q) 


2r'. 


The  elastic  equivalence  relations  were  then  written  as 


(13) 


tg  =  ^\ld 


(14) 


3Er(1-cos  (j)Q)w 

2^^3il  L 


The  general  buckling  stability  of  the  radome  was  evaluated  by  applying 

[21 

the  Tsien-VonKarman  equation ■*  to  the  equivalent  spherical  shell.  Substitu¬ 
tion  of  the  equivalence  relations  into  the  buckling  equation  defines  a  critical 
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t 


buckling  pressure  in  terms  of  the  element  dimensions  which  is 

.2 


(15) 


P  = 
cr 


3KEj^|i(l-cos 
2^/3i^I  R^L 


The  axial  stress  developed  in  the  space  frame  elements  was  evaluated 
by  applying  the  general  membrane  theory  of  shells.  The  partial  differential 
equations  of  equilibrium  for  the  equivalent  shell,  subjected  to  a  wind  pressure 
distribution  are 


(16) 


(.in  i  N^)  t  R^(N^g)  -  RNgCO.  ^  =  0 


(17) 


a 


R5^(.m  t  +  RN^g  cos  ^  =  0 


(18) 


Ni  +  N-  =  -  P„R  sin  4  cos  0 
f  0  0  ^ 


where  the  angles  ^  and  9  are  defined  in  Figure  2.  These  equations  have 

fsl 

the  general  solutions^-* 


(19) 


cos  0 

•  3* 
sin  0 

2 

•cos  ^  +  P  R(co8  ^  -  •:^os  ^) 
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'  5-  - 


„  sin  0  ^1 

♦»  "  ~ 

'  Sin  f 


'2  ^2 

—  +  - ^ - cos  ^  +  PqR(cos  ^ 


^  3x, 

-  ^os  f ) 


Nb  =  -(n^  +  p  qR  sin  ^  cos  ^ 


In  order  to  obtain  unique  solutions  for  the  stress  resultants  (N^,  N^,  N^g) 
equations  (19)  and  (20)  were  made  to  satisfy  two  boundary  conditions  and  the  in¬ 
tegration  constants  and  were  determined. 

The  first  boundary  condition  was  that  the  eicternal  mor/B  nt  about  an  axis 
in  the  plane  of  the  radome  base  and  normal  to  the  wind  vector  must  be  equal  to 
the  moment  due  to  the  internal  reaction  forces  of  the  shell  about  the  same  axis. 
Since  the  axis  is  in  the  plane  of  the  base  the  only  stress  resultant  contributing 
to  the  reaction  moment  is  N^. 

The  external  moment  is  caused  by  the  x  component  of  the  normal  wind 
pressure  (see  Figure  2)  and  is  equal  to 

2ir  ^0 


=  I  ^  PgR^cos  ^gSin^^cos^O  dj>  dO  . 


0  0- 


The  reaction  moment  due  to  the  meridian  stress  resultant  at  the  base 


(Ni  )  equals 


M.  =  \  Ni  cos  S  d9 

‘  ’^'>1 


where  Nx  indicates  the  stress  resultant  evaliiated  at  the  radome  base. 

"b 

Integrating  and  equating  the  two  expressions  gives  the  following  equation 
which  sat'.sfi.'^s  the  first  boiuidary  condition 


(24) 


7^  =  iPjR  cos  -  CO.  ♦jsm\  +  2)3  . 


The  second  boiuidary  condition  was  that  the  total  horizontal  external 
force  be  equal  to  the  total  horizontal  reaction  at  the  radome  base.  Two  stress 
resultants  contribute  to  the  horiz'^ntal  reaction,  and  The  total 

horizontal  external  force  equals 

Ztt  ^b 

r 


(25) 


t  007 

=  I  RpR  sin  ^  cos  d  d^  dO. 


xe  1  i  0 
0  0 

The  total  horizontal  reaction  at  the  base  equals 

Ztt 


2ir 


(26) 


F^i  =  R  sin  j  sin  0  dO  +  R  sin^^^cos  cos  0 


dO 


where  the  stress  resultants  and  are  evaluated  at  the  radome  base. 

Integrating  and  equating'  the  two  expressions  gives  the  following  equation 
which  satisfies  the  second  boundary  condition. 


(27) 


^6  0 

cos*  . 


cos  0 


sin  0 


sin 
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The  simtiltaneouB  solution  of  equations  (24)  and  (27)  provide  the  solu¬ 
tion  for  the  integration  constants  as  a  function  of  the  meridian  base  angle  A  . 
Qne-half  the  sum  and  difference  of  the  integration  constants  are  plotted  in 
Figvire  3. 

A  secondary  axial  stress  is' developed  in  the  equivalent  shell  due  to  the 
weight  of  the  structure.  The  equilibrium  equations  for  a  spherical  shell  sub¬ 
ject  to  a  imiform  gravity  load  are 


(28) 


R^(8iA  ^  N^)  +  (N^g)  -  RN^cos  0  =  tgCOS  ^ 


8 


(29)  ^  ^  RN^gCOS  e  =  0 


8 


(30) 


=  px.  sin  ^cos  0 


The  solution  of  these  equations  for  the  stress  resultants  are  given 


as 


W 


(31) 
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(33) 
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The  total  axial  stress  resultants  are  equal  to  the  combination  of 
each  component  due  to  the  external  loads  of  wind  pressure  and  structural 
weight. 

The  axial  stress  that  the  space  frame  elements  must  sustain  is 
related  to  the  stress  resultants  produced  in  the  equivalent  shell.  The  maxi* 
mum  element  stress  consists  of  the  product  of  the  effective  spacing  of  the 
elements  and  the  maximum  compressive  stress  resultant  divided  by  the 
cross  sectional  area  of  the  element,  or 


(34) 


2^^ 3(N  )L 

_ max 

^^'^ax  3(1 -cos  ^j^)wd  * 


The  maximum  compressive  stress. resultant  exists  at  the  wind  stag¬ 
nation  point  where  0  =  0  and  ^  =it/2.  The  general  stress  resultant  equations 
simplify  to 


(35) 


Na  =|(C,+ CJ  -^Rt„ 
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(36) 
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(37) 


at  the  stagnation  point  and  become  the  principle  axial  stress  resultants. 

The  v/ind  pressure  distribution  on  the  surface  of  the  radome  transmits 
a  transverse  force  to  the  individual  elements.  This  force  distribution  on  the 
elements  causes  flexure  stress  to  be  developed  in  addition  to  the  axial  stress 
expressed  in  equation(34).  In  analyzing  the  flexure  stress,  the  transverse 
force  distribution  was  assumed  triangular  with  a  maximum  unit  force  located 
at  mid-span.  A  distribution  was  assumed  in  order  to  eliminate  contention 
with  the  elasticity  of  the  dielectric  panels  and  the  consideration  of  elastic 
foundation  theory. 

The  magnitude  of  the  transverse  force  distribution  was  d^ned  such 
that  the  total  maximum  force  resisted  by  the  element  eqxialled  two-thirds  of 
the  total  external  force  applied  to  the  entire  panel.  The  distribution  varies 
from  zero  at  the  ends  of  the  element  to  a  maximum  at  mid-span  and  was  ex¬ 
pressed  as 
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^ 3P„* , 


o<x<i 


(38) 


q  = 


-  •/'S  Pq(x  -~), 
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The  maximum  flexure  stress  developed  from  simple  beam  theory  for 
a  condition  element  end  fixity  equalst 
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(39) 


,  6M 

(T^  -  max 

f  max 


wd 


16  ^0  ,2 
wd 


The  design  criterion  adopted  which  governs  the  selection  of  element 
dimensions  when  combination  of  column  (axial  compression)  and  beam  (flexure) 
loading  occurs,  was  such  that  the  stun  of  the  ratios  of  actual  to  allowable  axial 

The  allowable  flexure 
stress  at  the  outer  fiber  of  the  element  was  chosen  to  be  the  material  yield 
stress  and  the  allowable  axial  stress  was  defined  as  the  critical  Eulu  buckling 
stress.  The  design  criterion  was  expressed  as 


fiol 

and  flexure  stress  be  equal  to  or  less  than  unity. ' - ' 


(40) 


(T  r 

f  max  ^  max 


yp 


<  1 


which  upon  the  substitution  of  the  maximum  flexure  stress  equation  (39)  and 
Eulu's  critical  stress  equation  (40)  becomes 


(41) 
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The  development  of  equations  (15),  (39),  and  (41)  completes  the  deri¬ 
vation  of  the  structural  design  constraints.  The  constraint  equations  were 
developed  within  the  framework  of  simple  beam  theory  and  the  shell  membrane 
theory.  It  is  shown  in  ref[3]  that  the  ratio  of  membrane  stress  to  bending 
stress  is  small.  Therefore,  the  assumption  that  the  radome  equivalent  shell 
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behaves  like  a  membrane  is  valid,  provided  that  the  base  edge  of  the  radome 
is  free  to  expand.  However  the  base  of  the  radome  is  fixed  by  a  restraining 
ring  causing  a  discrepancy  in  the  boundary  conditions.  This  situation  is  recti¬ 
fied  by  superimposing  a  distributed  moment  around  the  edge  of  the  shell  of  such 
a  magnitude  that  the  membrane  displacements  are  made  to  vanish.  The  addi¬ 
tion  of  a  reaction  moment  to  the  base  results  in  an  increase  of  stress  in  the 
neighborhood  of  the  base  and  element  reinforcement  design  is  required  to  ac¬ 
commodate  this  condition. 

The  calculation  of  the  reaction  moments  at  the  radome  base  requires 
consideration  of  the  general  theory  of  bending  in  shells.  A  formulation  and 
solution  of  this  problem  is  available  in  refj^Sj ,  section  130,  which  gives  two 
particular  equations 


(42) 
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where  and  are  the  displacement  and  rotation  of  the  radome  edge 

and  Ml  and  H  are  the  reaction  moment  and  force  at  the  radome  edge. 

”o 


(44) 
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In  order  to  utilize  these  equations  for  the  solution  of  the  reactive  moment 
at  the  base  the  displacement  at  the  base  must  be  evaluated. 

The  computation  for  the  displacement  is  greatly  simplified  if  the  external 
loading  is  assumed  symmetrical.  Therefore  a  modified  external  load  is  imposed 
on  the  radome  to  facilitate  the  calculation  of  the  reaction  moment.  The  loading 
is  assumed  to  be  a  combination  of  a  uniform  normal  pressure  equal  in  magnitude 
to  the  wind  stagnation  pressure  and  a  uniform  structural  weight. 

The  displacer-'icnt  is  shown  in  the  refjsj  to  be  equal  to 


(45) 
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f  where 
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By  substitution  of  equations(46)  and  (47)  into  equation  (45)  the  displacement 
becomes 


(48) 


R  sin 
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The  boundary  conditions  of  radome  base  consist  of  reqmring  the  rotation 

to  equal  zero  and  the  displacement  to  be  equal  and  opposite  of  .  Therefore 

"b 

by  setting  equal  to  zero  and  substituting  equation  (48)  into  equation  (42) 

the  reaction  at  the  radome  base  was  evaluated  directly  and  equals 


Rt, 


(49) 
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The  optimum  elements  which  satisfy  the  structural  constraint  expressed 
by  eqviation(41 )  must  be  reinforced  to  accommodate  the  reaction  moments  dis¬ 
tributed  around  the  base.  The  degree  of  reinforcement  was  determined  by  re¬ 
quiring  that  the  maximum  stress  in  the  modified  elements  be  no  greater  than 
the  maximum  stress  in  the  optimum  element. 

The  reinforcement  design  is  completed  by  specifying  the  number  of 
elements  to  be  modified.  It  was  necessary,  then,  to  investigate  the  nature  of 

r4i 

the  bending  solutions  in  the  neighborhood  of  the  base.  It  has  been  shown*-  -*  that 
the  bending  moment  equations  can  be  expressed  in  the  form 


(50) 
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(52) 


^4  2. 

fb  =  3(1  -  H  ) — Y'  • 

The  exponential  term  in  the  moment  equation  causes  the  rapid  decay 
of  the  moment  at  the  radome  base  as  the  angle  a  (see  Figure  2)  increases,  and 
by  evaluating  the  rate  of  decay  the  number  of  panel  tiers  requiring  reinforced 
elements  was  determined. 

This  completes  the  description  of  the  structural  analysis  which  provides 
the  equations  and  techniques  to  evaluate,  (1)  a  structural  constraint  for  the 
majority  of  elements  distal  to  the  radome  base,  (2)  a  procedure  to  calculate 
the  reinforcement  reqiiired  by  the  elements  proximal  to  the  radome  base,  and 
(3)  the  member  of  elements  requiring  reinforcement. 

OPTIMUM  DESIGN  CRITERION 

The  concept  of  optimum  design  implies  a  methodilogical  assignment  of 
system  parameters  such  that  the  system  performance  is  maximized  relative 
to  a  defined  criterion.  There  are  several  possible  criteria  which  could  be  ap¬ 
plied  to  radome  design.  Two  criteria  were  selected  for  consideration  in  this 
study  and  are  (1)  minimization  of  the  radome  transmission  loss  subject  to  a 
holonomic  structural  constraint,  and  (2)  the  reduction  of  the  total  element  weight 
subject  to  an  electrical  constraint  which  fixes  a  minimum  transmission  loss  and 
a  structural  constraint  insuring  radome  stability  and  strength. 

A  description  of  the  optimization  procedures  initially  requires  a  defini¬ 
tion  of  a  single  structural  constraint  based  on  the  equations  developed  in  the 
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structural  analysis.  The  structural  constraint  utilised  was  the  most  severe 
design  condition  imposed  on  the  element  sizing  by  either  the  general  or  local 
buckling  equations  expressed  by  equations  (15)  and  (41).  A  comparison  of  the 
two  design  equations  was  made  by  graphing  the  maximum  allowable  wind  velocity 
permitted  by  each  equation  (see  Figures  4,  5,  and  6).  The  equation  which  per> 
mits  the  least  wind  velocity  for  a  particular  set  structural  and  environmental 
conditions  becomes  the  constraint  and  is  written  in  the  form 


F(L,  w,  d)  =  0 


The  equation  for  the  transmission  loss  is  also  a  function  of  the  element  dimen¬ 
sions  and  is  written  as 


G  =  G(L.  w,  d). 


The  procedure  for  obtaining  the  element  dimensions  which  will  result 
in  an  extremum  for  the  transmission  loss  (G)  subject  to  the  constraint  (I)  in¬ 
volves  the  application  of  the  method  of  Lagrangian  Multipliers  and  is  out¬ 

lined  below. 


The  total  differential  of  the  transmission  loss  is  written  and  set  equal 


to  zero: 


9G  9G  9G 
dG  =  g-^dL  +  dw  +-g^d  =  0 
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The  total  differential  of  the  structural  constraint  is  written 


(56) 


8F  3F  9F 
dF  =|^dL  +  -|^dw  +  r5dd  =  0  . 
dL  dw  dd 
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The  differentials  of  L,  w,  and  d  can  not  all  be  independent.  Arbitrarily 
specifying  the  differentials  of  L  and  w  as  indep^dent  the  remaining  dif¬ 
ferential  (dd)  is  solved  for  from  eqviation  (56)  by  setting 


(57) 


dd=  - 


|£dLt|:£dw) 

dF 
d  d 


provided  9F/9d  is.  not  equal  to  zero.  By  substituting  equation  (57)  into 


equation  (55)  the  dependent  differential  is  eliminated  and  the  resulting  equation 


is 


(58) 
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where  for  brevity  9G/9d)  /(9F/9d)  is  defined  as  Equation  (58)  can  be 
written  by  collecting  like  terms  as 


(59) 


9G  9f  ^F. 


Since  the  differentials  of  L  and  w  were  defined  as  independent,  the  coeffi¬ 
cients  must' vanish  in  order  to  satisfy  equation  (59).  In  this  manner  a  system 
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of  four  independent  equations  in  four  unknowns  (L,  w,  d,  \)  are  obtained,  the 
simultaneous  solution  of  which  will  yield  the  element  dimensions  that  cause 
the  transmission  loss  to  be  an  extremum.  The  transmission  loss  must  be 
evaluated  for  each  unique  set  of  element  dimensions  to  determine  if  a  mini- 
mxim  exists.  The  system  of  equations  to  be  solved  is 


(60) 


9<G 

8L 
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(61) 


(62) 


a.G 

8d 


=  0 


(63) 


F(L,  w,  d)  =  0  . 


A  direct  extension  of  this  procedure  would  consist  of  imposing  an 
additional  constraint  on  the  method  of  Lagrangian  Multipliers.  The  additional 
constraint  would  conceivably  be  placed  on  the  maximum  weight  of  the  space 
frame.  In  such  a  case  five  independent  equations  in  five  unknowns,  the  three 
element  dimensions  and  two  Lagrangian  Multipliers  and  would  be  solved 
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simultaneously.  A  complete  description  of  the  method  of  L»agrangian  Multi¬ 
pliers  can  be  found  in  ref[ll3  or  most  any  text  on  advanced  calculus. 

An_alternative  optimization  procedure  based  on  the  reduction  of  the 
structurtil  weight  subject  to  an  electrical  and  structural  constraint  was  de¬ 
veloped  and  is  essentially  graphical  in  complexion.  The  method  description 
is  first  precluded  by  the  formulation  of  the  required  equations.  The  structural 
constraint  will  be  the  severest  design  condition  imposed  by  either  the  general 
or  local  buckling  equations.  The  electrical  constraint  is  established  by  substi¬ 
tuting  a  value  fiar  the  minimum  allowable  transmission  loss  into  the  gain  equa¬ 
tion,  G  =  G(L,  w,  d).  The  electrical  constraint  equation  is  plotted  in  Figure  7 
for  constant  lengths  and  a  maximum  transmission  loss  of  0.  5db.  The  remaining 
equation  to  derive  is  the  expression  for  the  total  element  weight  in  terms  of 
element  dimensions.  This  equation  was  written  by  taking  the  product  of  the 
total  volume  and  the  weight  density  of  the  structural  elements.  The  total  vol¬ 
ume  is  equal  to  the  element  cross  sectional  area  times  the  collective  length 
of  the  elements.  The  weight  equation  was,  then,  expressed  as: 

N 

(64)  W  =  /^wd^L. 

} 

where  the  summation  of  element  lengths  has  been  shown  to  be  equal  to 
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Substituting  equation  (iO)  into  equation  (65)  gives 


(65) 


W  =  l3.85irR^ky>-~  . 


A  graphical  exanriiiation  of  the  three  equations,  electrical  and  structural 
constraints  and  the  structural  weight,  provides,  perhaps  a  less  direct  approach 
to  the  optimum  selection  of  element  dimensions.  However,  the  necessary  solu¬ 
tion  of  a  system  of  equations,  which  may  be  tedious  and  cumbersome,  is  avoided. 

The  objective  of  the  graphical  approach  was  to  generate  a  set  of  curves 
which  display  the  behavior  of  the  v/eight  and  constraint  equations  for  a  practical 
range  of  element  dimensions  and  arrive  at  a  design  decision  by  graphical  exami¬ 
nation.  The  curves  generated  for  the  example  design  were:  (I)  the  maximum 
wind  velocity  permitted  by  the  local  buckling  equation,  (2)  the  maximum  wind 
velocity  permitted  by  the  general  buckling  equation,  (3)  a  comparison  of  local  and 
general  buckling,  (4)  the  electrical  constraint  equation  for  constant  lengths, 

(5)  the  weight  of  aluminum  elements  for  a  range  of  widths  which  satisfy  both  an 
electrical  and  structural  constraint,  (6)  the  weight  of  aluminum  elements  for  a 
range  of  lengths  which  satisfy  both  an  electrical  and  structural  constraint,  (7) 
the  exponential  decay  function  for  the  bending  moments  due  to  the  radome  bound¬ 
ary  conditions. 
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SAMPLE  RADOME  DESIGN 

An  example  design  of  a  150  foot  diameter  space  frame  radome  was 
performed  to  provide  an  illustration  and  detailed  description  of  the  optimi¬ 
zation  design  procedures.  The  design  calculations  were  based  on  the  follow¬ 
ing  radome  and  environmental  specifications: 

1.  Radome  diameter  -  150  feet. 

2.  Structural  description  -  triangulated  space  frame  of  regular  (snub-dodeca¬ 
hedral)  or  random  (icosahedral)  geometry. 

3.  Truncation  base  diameter  -  130  feet. 

4.  Structural  material  -  6061  T6  aluminum. 

5.  Maximum  wind  velocity  -  150  mph  sustained. 

6.  RF  band  -  200-500  me. 

7.  Maximum  transmission  loss  -  0.  SOdfa. 

The  general  radome  design  equations  developed  in  the  preceding  sections 

are: 

1.  Transmission  loss  -  gain  equation 

C  -  ^  Tv/d  18d 

^  ^  2Lw 

2.  General  buckling,  Tsien-VonHarman  equation 
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3KE^(1-co8  t  )  wd^ 

~  - 

i 

3.  Local  buckling  criterion  ^ 

I 

> 

^^3P.L^  2^/3(P  R-VV/  Rd)L^  j 

- i! -  +  - 9 - ? -  <  i  I 

16(T'  wd^  ir^E(l-cos  4.  )wd^  ? 

yp  '  Tb  , 

4.  Total  weight  of  elements 

W  =  13.85irR  k/-^ 

These  equations  become  specific  design  formula#  upon  substitution  of 
the  appropriate  parameters  which  are: 

1.  =  10^  psi 

2.  K  =  0.50 

3.  k  =  0.75 

4.  ^Ij^  =  120° 

5.  R  =  900  inches 

6.  ^  =  40i  000  psi 

yp 

7.  /  =  0.  098  pcf 

8.  Q  =  0.  004  pcf 

9.  Pq  =  0.3125  psi 

10.  G  =  .  5db 


101 


The  design  formulae  are  equal  to 


(66) 


Lw  =  79  w 


2 


+  7wd  18d 


(67) 


P 

cr 


4.  64 


wd 

L 


(63) 


(1.316  +  0.  0832d)L^  x  10“^  =  wd^ 


(69)  W  =  2643^  (kips) 

The  sample  design  description  includes  the  method  of  Lagrangian 
Multipliers  to  minimize  the  transmission  loss  subject  to  one  constraint  and 
a  graphical  method  to  reduce  the  total  weight  of  elements  subject  to  an  elec¬ 
trical  and  structural  constraint. 

In  both  illustrations  it  is  necessary  to  eliminate  one  of  the  structural 
design  formulae  by  determining  which  dictates  the  most  severe  restriction  on 
the  element  dimensions.  Figures  4  and  5  are  plots  of  equations(67)  and  (68). 
Figure  6  is  a  comparative  plot  of  equations  (67)  and  (68)  and  displays  the  fact 
that  for  all  element  widths  and  lengths  over  one  inch  and  14  feet,  recpcctively, 
the  local  buckling  design  restriction  is  most  severe.  Since  a  width  of  one  inch 
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or  less  is  impractical  for  the  size  radome  considered  in  the  example  design, 
the  local  buckling  equation  was  selected  as  the  structural  constraint. 

The  method  of  minimizing  the  transmission  loss  subject  to  a  single  con¬ 
straint  was  previously  outlined  where  equations  (53)  and  (54)  take  on  the  form: 

(53)  F  =  (1.  316  +  0.  0832d)L\  lO”^-  wd^=  0 


(54) 


G  = 


79w  +  7dw  +  18d 
2L.W 


The  system  of  equations  to  be  solved  for  the  element  dimensions  which 
minimize  the  transmission  loss  are  obtained  by  performing  the  appropriate 
differentiation.  The  equations  in  terms  of  the  four  unknowns  L,  w,  d,  and  X 
are: 


(55) 

(56) 

(57) 

(58) 


39.  5w^+  3.  5dw  +  9d  +  3Xwfl.  316+  0.  0832^  L^x  10'^=  0 
39.  5w^-  9d  +  Xw^d^L  =  0 
3,  5w  +  9  -  Xw[0.  0832L^-  2wd^jL  x  J0"^=  0 
[l.  316  +  0.0832djL\  iO"^-  wd^=  0  . 


103 


The  method  of  Lagrangian  Multipliers,  which  requires  the  simultaneous 
solution  of  the  above  equations,  is  a  direct  analytic  approach  to  the  evaluation  of 
the  element  dimensions  consistent  with  minimum  transmission  loss.  The  above 
system  of  equations  was  developed  as  an  illustration  of  technique  only,  and  the 
solution  was  not  attempted.  However,  attention  is  directed  to  the  degree  of  al¬ 
gebraic  complication  involved  in  obtaining  a  complete  solution  which  would  be 
increased  by  the  addition  of  a  second  constraint.  Therefore,  emphasis  was  cen¬ 
tered  on  a  graphical  approach  which  although  less  analytically  direct,  is  more 
favorable  to  hand  calculation  and  visual  interpretation. 

The  graphical  method  w.as  instrumented  by  the  generation  of  a  set  of 
curves,  Figures  4  through  9.  The  total  weight  versus  element  width  curve, 
Figure  8.  clearly  defines  a  consequential  design  objective  which  parallels  the 
reduction  of  weight,  which  is  the  minimization  of  width.  However,  it  is  neces¬ 
sary  to  establish  a  minimum  width  below  which  the  advantage  of  suppressed 
structural  weight  is  totally  masked  by  vinmanageable  mechanical  design  difficul¬ 
ties.  The  total  weight  versus  element  length  curve.  Figure  9.  provides  a  cue 
to  the  selection  of  a  minimum  element  width.  The  weight  equation  is  not  a  single 
valued  function  and  displays  a  point  of  minimum  element  length  which  satisfies 
the  electrical  and  structural  constraints.  The  element  width  which  corresponds 
to  the  minimum  element  length  is  of  a  practical  magnitude  and  is  defined  as  the 
goal  for  the  ensuing  mechanical  design.  Therefore,  the  requirements  for  struc¬ 
tural  integrity  during  sustained  150  mph  winds,  minimum  transmission  loss  of 
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0.  50db>  minimum  element  length,  and  minimum  total  weight  of  elements  within 
practical  limitations  are  satisfied  for  a  structural  element  of  the  following  di¬ 
mensions: 

w  =  1. 10  inch 
d  =  3.42  inches 
L  =  167  inches 

The  remaining  design  consideration  concerns  the  reinforcement  of  the 
standard  element.  The  reinforcement  calculation  involved  the  evaluation  of  the 
shell  displacement  and  derivatives  at  the  radome  base  and  substitution  into 
the  bending  moment  equations  (48)  and  (49)  which  gives  for  the  resultant  mo¬ 
ment  at  the  base: 


(70) 


^0 


PoRd 


3(l-tii) 

p.(l-|i 


The  bending  moment  applied  to  the  elements  at  the  base  are  then  eqxuil  to  the 
resulting  moment  Mx  times  the  effective  spacing  (s). 

Tfl 


(71) 


Mx  s  =  i5.95dL. 
^0 


The  standard  element  was  modified  such  that  the  maximum  flexure 
drcBB  of  the  reinforced  elements  was  maintained  equal  to  the  maximum  flexure 
stress  of  the  standard  element.  This  reqxiirement  was  formulated  by  equating 
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the  flexure  stress  of  the  two  elements,  one  of  standard  cross  section  and  the 
other  of  modified  cross  section. 


(72) 


'vTs  P  L" 
0 

16wd^ 


V- 3 


6M.  s 


16(w+/^)(di-  {w+^w)(di-^d)^ 


Equation  (72)  was  solved  explicitly  for  the  incremental  width  ^w  and  equals 


(73) 


wd^  +  (1530-1-  90d)wd^ 

(cH-^d)^  3  P^L^(df^d)^ 


The  most  convenient  manner  to  reinforce  the  element  is  to  maintain  a 
fixed  element  depth  and  increase  only  the  width.  By  fixing  the  element  depth 
for  the  entire  space  frame  the  requirement  for  special  transition  hubs  or  adap¬ 
ters  is  avoided.  The  reinforcement  equation  is  simplified  by  letting  ^d 
vanish  and  reduces  to! 


(74) 


^w  =  (2830-H  167d) 


wd 


The  requirement  in  element  width  for  the  sample  design  equals  0.  459 
inch  and  the  reinforced  element  dimensions  become 
w  =  1 . 559  inch 
d  3.  42  inch 
L  =  167.  0  inches 
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The  number  of  elements  requiring  reinforcement  was  determined  by 
examining  the  rate  at  which  the  bending  moment  at  the  radome  base  damps 
out  at  points  away  from  the  base.  The  percent  reduction  in  the  moment  was 
plotted  as  a  function  of  meridian  angle  (a)  receding  from  tl;e  base  as  dictated 
by  the  exponential  term  of  equation  (50).  Figure  10  shows  that  for  the  sample 
design  the  moments  are  reduced  by  98%  a.t  an  angular  displacement  of  10°  from 
the  base.  An  equation  for  the  number  of  reinforced  elements  was  derived  in 
the  following  manner: 

Let  A^  be  the  total  area  of  the  radome  composed  of  reinforced  ele¬ 
ments.  Then  A^  is  eqvial  tot 

(75)  A^  =  ZtrR^a  sin{^^-  |)  . 


The  number  of  reinforced  panels  (M^)  equals  the  area  divided  by  the  unit 
panel  area  and  equals 


(76) 


M 


8trR  a  sin(^  ~  ~) 

U  Cm 


^/■3  V 


and  the  number  of  elements  corresponding  to 


M 

r 


panels  is 


(77) 
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or 


(78) 


N  = 
r 


12irR  a  sin(^j^- 


The  value  of  for  the  sample  design  and  a  98%  reduction  of  shell 
induced  bending  becomes  66  elements.  The  total  number  of  elements  of  both 
cross  sections  is  determined  by  substituting  the  element  length  into  equation 
(11),  which  gives 


N  =  977. 


The  complete  specification  of  the  elements  of  the  sample  design  is 

1.  Standard  Element 

Number  -  911 
Width  -  1, 10  inches 
Depth  -  3. 42  inches 
Length  -  167,  0  inches 
Weight  -  55,400  lbs 

2.  Reinforced  Elements 

Number  -  66 
Width  -  1.48  inches 
Depth  -  3.42  inches 
Length  -  167.0  inches 
Weight  -  5,  594  lbs 
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MECHANICAL  DESIGN  OBJECTIVES 


There  are  numerous  mechanical  design  deficiencies  existing  in  many 
of  the  large  radomes  which  have  been  built  in  the  recent  past.  This  fact  was 
made  apparent  during  consultation  with  the  personnel  of  several  prominent 
manufacturers  cf  large  space  frame  radomes  as  well  as  by  the  inspection  of 
the  Haystack  Installation.  The  most  blatant  deficiency  seems  to  be  in  the  area 
of  dielectric  panel  to  element  attachment.  A  result  of  manufacturing  panels 
of  reinforced  plastics  within  reasonable  tolerances  seems  to  be  a  condition  of 
relative  slack  in  the  individual  panels'  membranes  after  assembly.  Such  a 
condition  allows  considerable  furling  of  the  panels  due  to  wind  pressure  fluctu¬ 
ation.  Fatigue  failures  of  the  panels  frequently  occur  at  the  panel  vertices 
which  require  undue  and  excessive  maintenance.  A  possible  solution  or  sup¬ 
pression  of  this  type  problem  wovild  be  effected  if  an  attachment  design  was 
generated  which  produced  a  prehension  in  the  panel  during  installation.  Such 
a  design  objective  is  complicated  by  the  optimum  design  conclusion  that  mini¬ 
mum  element  width  is  desirable.  A  design  concept  is  depicted  in  Figure  11 
which  has  the  merit  of  simplicity  and  warrants  development  consideration.  The 
extruded  V-slot  in  the  edge  of  the  element  would  cam  the  reinforced  panel  taut 
during  depression  of  the  panel  into  the  slot  by  the  cover  plate.  The  camming 
action  would  induce  a  prestress  on  the  panel  dependent  on  the  pitch  of  the  slot 
and  depth  of  depression. 
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A  second  deficiency  seems  to  exist  in  obtaining  and  maintaining  ade¬ 
quate  weather  seals.  The  panel  to  element  seal  coxild  be  significantly  improved 
if  and  when  prestressed  panels  become  reality  due  to  the  increased  bearhig 
pressure  and  area  of  contact  at  the  seal.  The  panel  to  hub  seal  is  generally  the 
most  difficult  to  obtain.  By  providing  a  panel-hub  engagement  slot  in  the  hub 
a  forced  and  continuous  attachment  of  the  panel  to  the  element  and  hub  is  possible. 
Therefore,  direct  contact  of  a  hub  cover  plate  and  the  panel  would  allow  for 
generous  surface  for  a  continuous  weather  seal.  Figure  12  schematically  il¬ 
lustrates  such  a  panel-element-hub  attachment  scheme. 

The  preceding  recommendations  are  intended  only  as  conceptual  pro¬ 
posals  which  could  be  evaluated  for  def»ign  feasibility  and  refinement.  A  multi¬ 
tude  of  design  improvements  in  the  mechanical  design  of  radomes  all  of  which 
require  considerable  analysis.  Some  of  the  areas  for  improvement  include 
erection  procedures,  reduction  of  manufacturing  tolerances  to  a  minimum  for 
mating  radome  components,  elimination  of  superfluous  assembly  hardware, 
and  the  aforementioned  panel  pretension  and  weather  seals  each  of  which  re¬ 
quires  individual  and  special  attention. 
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APPENDIX 


DERIVATION  OF  EXPECTED  SIDE  LOBE  LEVEL  DUE  TO  RADOME 

From  equation  (24)  of  QJ  we  may  write  the  far  field  pattern  in  the 
presence  of  a  space  frame  as 

Zir  a 

(Al)  F(u)  =  I  I  r  dr  (1  + g^^^(r,/|  ))f(r, 

0  0 

where  f(r,^)  is  the  antenna  aperture  illumination,  is  the  blocking 

area  ratio  considered  as  a  function  of  the  polar  coordinates  r,  (see  p.  12 
of  £lj )  and  u  =  k  sin  0  where  0  is  the  pattern  observation  angle.  We  shall 
now  find  the  expected  value  of  the  power  pattern  jF(u)j  under  the  assumption 
that  is  one  of  a  large  set  of  functions  denumerated  by  a  variable  a 

of  known  statistics  approximating  those  of  typical  space  frames.  Hence  we 
write  as  y7(r,/|,a).  The  average  value  of  y^*(r,/^,a)  is  the  average 

blocking  area  ratio  ^  so  that 

(A2)  ( y?(r,  f|,a)da  =  p,  where  fda  =  l 


and  the  integrals  are  over  all  admissible  a,  normalized  according  to  (A2). 

We  may  then  write  the  expected  value  of  the  power  pattern  from  (Al) 
as 
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(A3) 

a  2ir  a  2ir 

+  jpdp  I  drj  |r  dr  |d^  1  f(r,r|^)£(p,  cos  -  p  cos 

0  0  0  0 

where  ^Q{’i)  is  the  unperturbed  pattern  and  ihe  inner  integral 
(A4)  I--  |^day)(r,r^,a)y>  (p,  ^,a) 

may  be  computed  for  any  given  space  frame  statistics.  On  the  basis  of  some 
detailed  computations  a  good  approximation  is  believed  to  be  Gaussian  with  the 
scale  length,  the  width  of  the  element.  It  is  quite  obvious  that  two  points  more 
than  W  apart  are  not  highly  correlated. 

(A5)  I  ^  )e  -  /^e 

_  _  2  2  2 

where  [r-pj  =r+p-  2rp  cos  (♦-a)- 

If  we  make  the  further  reasonable  assumption  that  in  the  neighborhood 
of  appreciable  correlation,  |r  -  p|  <  ‘~W,  is  virtually  constant,  we 

may  explicitly  integrate  the  last  term  in  (A3)  to  obtain 
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t 


,  ,  ,  ,  ,  ,  -WW/Z)^ 

(A6)  jiF(u)|  da  =  I !  t  |Fu(u  )j  +  yW  |g^J  e 

Zt!  a 

I  2 

^  rf  (r,  ^)dr. 

0  0 

The  increase  of  pattern  in  the  neighborhood  of  the  first  few  side 
lobes,  where 

-{uW/2)^ 

(A7)  e  ^1. 


relative  to  the  peak  is  given  then  by 


(A8) 


2i 

I  jdi  ^  rf(r,  ^)drj^ 
0  0 


r 


The  expression  (Z8)  in  the  text  is  an  approximation  to  (A8)  using  (2)  and 
typical  values  of  f(r,  ^). 


1 


119 


DERIVATION  OF  EQUATION  (1)  Part  1 


The  far  field  of  an  aperture  anter.na  may  be  written  as  the 

Fourier  transform  of  its  aperture  illumination 


(A9) 


Fq  =  Fq{0)  =  (  [  f(r.  ®  r  dr  da 


when  0  is  a  unit  vector  in  the  direction  of  observation,  r  is  a  vector  from 
the  origin  to  the  integration  point,  k  is  the  wave  number,  and  r  and 
are  polar  coordinates  in  the  aperture. 

The  induced  field  ratio  IFR  of  a  set  of  scatterers  in  the  antenna 
beam  (1)  whose  mutual  coupling  can  be  neglected  and  (2)  whose  dimensions 
projected  into  the  antenna  beam  are  small  compared  to  variations  in 
f(r.  ^|)e''  and  are  non-overlapping,  can  be  defined  as  g{r./]),  equal  to 

either  the  IFR  of  the  scatterer  projecting  to  (r.  /|  )  or  equal  to  zero  if  no 
scatterer  projects  to  r,  .  In  this  case  the  scattered  far  field  may  be 
written  as 


(AlO) 


F^  =  F^(0)  =  ^rdrd^. 


K  all  of  the  scatterers  have  the  same  IFR,  g,  a  constant,  then  (A  10) 
becomes 


(All) 


=  /g  ||  d/i. 


no 


where  p  is  the  ratio  of  the  projected  area  occupied  by  the  scatterers  to  the 
total  area  of  the  beam,  or  the  blocking  area  ratio.  From  (A9)  and  (All)  we 
may  w.-ite  the  total  far  field  as 


(A  12) 


In  the  forward  direction  (2)  above  is  satisfied  for  all  practical 
space  frames  and  the  relative  power  loss  due  to  the  radome  is 


(A13) 


U  +/gl^  • 
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